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Abstract 

In this article we obtain some general explicit formulas to express the 
n-th term of a wide class of sequences of numbers defined by recursion 
in terms of the partitions of n. These formulas are here applied to 
the Bernoulli numbers, to the Euler numbers and to the inverse of the 
factorials. Our proofs are based only on elementary linear algebra and 
do not use any tools from calculus. 



1 Introduction 

Many sequences of numbers are defined through a recursive formula. For 
example, Lucas in [1] (Eq. (3), p. 238, and Eq. (1) p. 256) gives the 
following symbolic equations defining Bernoulli and Euler numbers: 

(B + l) n -B n ^ 0, Vn>l, Bq = 1, #i = -~ 
(E + l) n + (E — l) n ^ 0, Vn > 0, E = l 

These equations are symbolic because after the expansion the exponents must 
be degraded to indices. The resulting recurrence relations are the following: 



h=0 



£ (f) [1 + (-irl E h = 0, E = l 



(2) 
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(where in Eq. ([Tj) n has been substituted with n+ 1, and the definition of Bi 
is no longer necessary) . One may then write the following recursive formulas 
for B n and E n : 



n-l 

I 

Br 



1 fn + 1 \ 

— £( ft * = i (»> 

n-l / x 

£ --oE j [i + HM^, E = l (4) 

Another way to define these sequences of numbers is through a generating 
function. In this case the numbers of the sequence are related to the co- 
efficients of the Taylor series of some generating function. For example, a 
possible definition of the Bernoulli and Euler numbers is following (pQ, p. 
262): 

00 h 
X \-— \ „ X 



h=0 
h=0 

Definitions like those shown have the disadvantage that to calculate the n- 
th element of the sequence one has to calculate all the elements with index 
smaller than n. 

Formulas that allow to calculate the n-th element of the sequence only 
in terms of n, and not in terms of other elements of the sequence, are then 
useful and welcome. They are often called closed formulas, but we prefer to 
call them explicit formulas, in agreement with [2]. 

In this article we present in Theorems [T] and |5] explicit formulas that 
allow to calculate the n-th element of a sequence of numbers in terms of the 
partitions of the integer n. This is possible if the recursive definition of the 
sequence elements is of a certain form, not very restricted, as we will see. 
We then apply Theorem Q] to find explicit formulas for the inverse of the 
factorials, the Bernoulli numbers and the Euler numbers. Our results partly 
overlap with those obtained by Vella [2] and by Malenfant j3], but we think 
that our proofs are more elementary because they are done just using some 
basic linear algebra and some combinatorics. 

Before to state our theorems it is necessary to give some definitions con- 
cerning partitions, compositions and Toeplitz matrices. 

Given a positive integer n, a partition p of n is an unordered sequence of 
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Kp) positive integers rij that have sum n: 

Kp) 

p = {ni,n 2 ,...,ni(p)}, n = ^rij 

3=1 

The number l(p) is called the length of the partition p. 

A composition c is a partition in which all the n.,- are ordered, for example: 

nj>n j+1 , V? = 1, . . .,Z(c) - 1. 

For any composition c there are many partitions formed with the same 1(c) 
elements of c. Let //(c) be the number of different partitions formed by the 
same numbers forming the composition c. It is an easy exercise to verify that 
the number //(c) is given by the formula: 

//(c) = /(c)! / J] m c ( nj )\ , (5) 
/ rij e u(c) 

where m c (rij) is the multiplicity of in c, that is the number of times that 
the number rij appears in c, and the product at the denominator has to be 
done for all different numbers rij that appear in c, that is, for all numbers 
rij that appear in the union U(c) of the elements of c. The multiplicities 
m c (jij) that appear at the denominator of §5§ always form a partition of the 
length 1(c) that appears at the numerator. If p is a partition obtained from 
the composition c, one obviously has l(p) = 1(c). 

Let C(n) and V(n) be the set of all compositions of n and the set of all 
partitions of n, respectively 

An n-dimensional Toeplitz matrix A is a square matrix of order n that is 
constant along all diagonals, that is, such that: 

Aij = Vi, j = 1,2, ...,n- 1. 

An interesting way to define a Toeplitz matrix is through a function 
defined on the set Z of the integers: One just has to define the elements Aij 
of a square matrix A as functions of the differences i — j, that is, one just 
has to define: 

A id = f A (i- j)- 

The index A in the function symbol j a is just to point out that the function f A 
completely defines the Toeplitz matrix A. A lower triangular Toeplitz matrix 
is a Toeplitz matrix with vanishing elements above the main diagonal, that 
is, such that 

A id = 0, Vi<j. 
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A non-singular lower triangular Toeplitz matrix has necessarily a non-zero 
number in the elements of its main diagonal. If the lower triangular Toeplitz 
matrix A is defined trough a function Ja, then A is non-singular if and only 
if /a(0) / 0. 

To simplify the notation, we will put 

a k = f A (k), Mk = 0,1,2,.... (6) 

We are ready to state our main results. 

Theorem 1 Letn > be any given positive integer and f a a function defined 
in the set of non-negative integers such that /U(0) 7^ 0. Let be defined by 
Eq. If a sequence a , at, . . . , a n , satisfies the following condition: 

n 
h=0 

which is equivalent to the following recursive definition of the elements a n of 
the sequence: 

^ n— 1 

a n = V 

one may calculate the elements a n , Wn > 0, from the following explicit for- 
mulas involving the partitions of n: 

a. = m £ II 0» 



n-1 



(10) 



ceC(n) rij-ec 



Theorem 2 Let ao,a%, . . . ,a n be a sequence to which Theorem^ is applica- 
ble, and so Eqs. (TjJ) and [ffl\! are true, then also the following formulas are 
tru e : 

a, <-!)■ = «. £ II i-^ 1 ) (") 

peV{n) n^ep v u y 

a n (-ir = a £ Mc) ft (-^^) ( 12 ) 
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If a Q = 0, Eqs. (15])- (fT2"]) give a n = 0, Vn > 0, so, to avoid trivial 
results, one may suppose a ^ 0. Then, by dividing by a all Eqs. f[T|)- f[T2"j) . 
one sees that only the fractions a n /ao are actually present in Eqs. (T7|)- (fT2|) 
and depend on the sequence {ccq, a%, . . .} (that is on /a). This allows to 
consider only sequences of numbers {ao,ai, ■ ■ ■} normalized in such a way 
that a = 1. Moreover, the function f^, or, equivalently all the elements 
of the sequence {a , act, . . .}, can be multiplied by a constant factor without 
changing anything in Eqs. fP7|)- ffl~2]) . so we may suppose in all generality, that 
f A (0) = a = l. 

Formulas ([9]) and ffTOl) in Theorem [1] or ffTTj) and ffT2]) in Theorem [2] differ 
only for the use of partitions in place of compositions. In general, the use 
of compositions is preferable in calculations because it involves fewer terms. 
We state all subsequent results in terms of compositions only, but there exist 
similar formulas concerning partitions, too. 

By applying Theorem [1] we obtain the following results: 

Theorem 3 For any given integer n > 0, one has the following explicit for- 
mula to express the inverse of the factorial of n in terms of the compositions 
of n: 

^ - e f{c> n (-£) ■ (") 

// n is even, one also has the following explicit formula to express it in terms 
of the compositions of n of length greater than 1: 



c 6 C(n) Uj € c 

1(c) > 1 



// n is odd, one has the following formula in terms of the compositions of n 
of length greater than 1: 







E *>n(-^)- <«) 



c 6 C(n) llj £ c 

1(c) > 1 



Theorem 4 For any given integer n > 0, one has the following explicit 
formulas to calculate the Bernoulli numbers B n and B 2n in terms of the 
compositions of n: 

§ - E mc) n (-^) • m 
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(2n - 1) B 



2n 



(2n)! 



e ^ n 

cGC(n) 



n,' £ c 



(2n i + 2)!y 



(17) 



Theorem 5 For any a«i> en integer n > 0, one has the following explicit for- 
mulas to calculate the Euler numbers E n and E 2n in terms of the compositions 
ofn: 



E„ 



n 



e Me) n 

ceC(n) 



n 7 - G c 



2d 



e ^ n 



i + (-if; 



1 



(2n)! ^ 11 V (2n 7 -)! 7 ' 



(18) 
(19) 



Let's show the use of Formulas f|T5]) - ffT^]) with some examples. 

Example 1. Some formulas using C(4). One has: 

C(4) = {{4}, {3, 1}, {2, 2}, {2, 1, 1}, {1, 1, 1, 1}} , 
so from (fl6|) one finds: 

'1! -1 2! (-1) 2 2! (-1) 2 



E 4 = 4! 



1! 5! 1!1! 4! 2! 2! 3! 3! 



+ 



+ 



3! 



•I) 3 4! 
+ 



1!2! 3! 2! 2! 4! 2! 2! 2! 2! 
while from (I14p one finds: 

-l) 2 3! (-1) 3 



1 

30 



1 
2 



2! 



;-l) 2 + 2! 



+ 



4! 



:-i) 4 



1!1! 3!1! 2! 2! 2! l!2! 2! 1! 1! 4! 1! 1! 1! 1! 



1 

24 



Example 2. Some formulas using C(5). One has: 

C(5) = {{5}, {4, 1}, {3, 2}, {3, 1, 1}, {2, 2, 1}, {2, 1, 1, 1}, {1, 1, 1, 1, 1}} 
so from f)16p one finds: 



B* = 5! 



1! -1 t 2! (-1) 2 t 2! (-1) 2 + ^|_ (-1) 3 | 



1! 6! 1!1! 5! 2! 1! 1! 4! 3! l!2! 4! 2! 2! 

= 0, 



3! (-1) 3 4! 

+ 777T7 • 777777777 + 



-I) 4 + 5! (-1) 5 



2!1! 3! 3! 2! l!3! 3! 2! 2! 2! 5! 2! 2! 2! 2! 2! 
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while from ( I15p one finds: 
2! (-1) 2 



+ 



1!1! 4!1! 1!1! 
3! f-1) 3 4! 



(-1) 2 
3! 2! 

-I) 4 



3! 



1!2! 3! 1! 1! 



5! 

+ ^7 



2!1! 2!2!1! 1!3! 2! 1! 1! 1! 5! 1! 1! 1! 1! 1! 
and from ffT9l) one finds: 



E 



10 



10! 



1! -1 
il ' 10! 



1!1! 



3! {-If 



2!1! 4! 4! 2! 



+ 



1!1! 8! 2! 

4! (-1) 4 5! 
l73! ' 4! 2! 2! 2! + 5! ' 2! 2! 2! 2! 2! 



6! 4! 

(-1) 5 



_3!_ (-1) 3 
l\2\ 6! 2! 2! 

= -50521 



+ 



Results similar to those stated by Eq. (fTB"j) and Eq. ( 1T9"j) were also found 
by Vella and by Malenfant. Vella in [2] obtained a formula for B n (Formula 
(a) of his Theorem 11) using the formula of Faa di Bruno |4J for the n-th 
derivative of a composite generating function. Malenfant in [3] obtained a 
formula for B n (his Formula (36)) and a formula for E 2n (his Formula (39)) 
starting from the coefficients of the Taylor series of a generating function. 
Our method is in some points close to that one of Malenfant but, contrary to 
his method, it does not use any tool from calculus. We think that the results 
expressed by Theorems [TH3] and by Eqs. (iT7j) and (fl8"j) are new. In any case, 
our proofs are simple because they only use basic tools of linear algebra and 
some combinatorics. 



2 Proof of Theorem H 



Let n be any fixed positive integer. Let's first take then out of the sum in 
Eq. ^ the contribution corresponding to h = 0, so Eq. (J7|) may be written 
in the following way: 

n 

2j a n - h a h = -a n a , (20) 

h=l 

Let's introduce the n- dimensional column vectors v and w whose elements 
are 

v k = a k) and w k = —a k a , 
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respectively, where k is the vector index running from 1 to n, and consider 
an order n square matrix A whose elements are defined by 

Aj = { V I>: ' , (2i) 

J l if i < j ' v 1 

so that A is a non-singular n-dimensional lower triangular Toeplitz matrix. 

Equation (I2(jp becomes now a linear system for the unknowns Vk whose 
matrix expression is the following: 

Av = w . (22) 

Being A non-singular, one formally writes the solution to the system (I22p in 
the following way: 

v = A~ l w . 

The element v n is then easily obtained from the scalar product of the vector 
at the last row of A' 1 with the vector w: 

n 

v n = ^2{A- 1 ) n , h w h . (23) 

h=l 

The inverse matrix A^ 1 could easily be calculated using any elementary 
method of linear algebra. It is itself a non-singular lower triangular Toeplitz 
matrix. There are also plenty of books and papers describing many algebraic 
and analytical properties of Toeplitz matrices (see for example [5], or [3]). 
Instead to follow this road, in order to arrive at the result expressed by The- 
orem (U we will proceed in a more elementary way, thanks to the particularly 
simple forms of A and w. For convenience let's put iV = A" 1 and write 
explicitly the matrix A: 



A = 



( 


OiO 








. 


. 


\ 




Oil 







. 


. 






« 2 


ox 


oca 


. 


. 






a 3 




Oil 




. 














. 




V 




O-n-2 


«n-3 


ftn-4 • 


• «o 


/ 



Equation fT23l) requires the knowledge of the last row of the matrix N. To 
get that we start from the equation: 

n 

$n,k = N nt h Ah t k 
h=l 
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and use the fact that A hjk = if h < k. We then find: 



for k = n : 1 = N n , n A n>n ; 

for k = n-h : = Y2=n-h N n,k A k ,n-h, V7i = 1,2, 1. 

The equations just written form a triangular linear system of n equations in 
the n indeterminates N n;h , h — 1, . . . , n, that is easily solved by substitution. 
By taking into account that A k>k = «o 7^ 0, VA; = 1, 2, . . . , n, one finds: 

1 

N, 



n,n , i 



n 



N n , n - h = -— 2^ N n , k A k , n _ h , V/i = l,2,...,ra-l. (24) 



k=n—h+l 



Each element N n ^h of the last row of the matrix N is then expressed in 
terms of those with greater column index. Let's explicitly solve the system 
for the first few values of h, substituting also the matrix elements of A with 
the elements of the sequence {ckq, «i, . . .}. One so finds: 



N, 



n,n 



1 _ J_ 

A n;n ao 



AT _ 1 AT A 1 I Ctl 

^n,n-l — Ann ^n,n^n,n-l— aQ y Q() 

N n , n -2 = ~- A - 



— (N n n -i A n _i n _ 2 + N n n A n n ^ 2 ) — ^- ( — 

^n,n— 3 = A n {^n,n— 2 A n — 2,n— 3 ,n— 1 -^n— l,n— 3 ,n ^4n,ri— 3 

_ J_ [ _ a l _i_ ^2^1 _i_ ai«2 _ Q3 | 
a \ ay a ( 2 , a 2 , «0 / 



1_ / £j Q2Q^ 010201 I a 3 oi a i a 2 _i_ a i _i_ 0103 04 \ 

~$ og a 2 q ° / 



At this point one notes that, at least for h — 1, 2, 3, 4, . . ., N n>n _ h is obtained 
as the sum of all products like 

whose indices satisfy the conditions 

1 < k < h, rii + n 2 + . . . + n k = h, 

and each term is taken with the + sign if k is even, with the — sign if k is 
odd. The k positive numbers rii, . . . ,n k form then a partition of length k of 
the positive integer h. Using the notation presented in the Introduction, we 
are then led to write the formula: 
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N n , n - h = — J2 IT—' Vfa=l,2,...,n-1, (25) 
a *— ' - LJ - On 

p&V(h) n d e P 

where V(h) is the set of all partitions of the integer fa and l(p) is the length 
of the partition p. Let's verify by induction that Formula ( 12 5 p holds true for 
all values of h — 1,2, . . . ,n — 1. Supposing (j25|) true up to a given fa, with 
1 < fa < n — 2, let's prove that it holds true also for fa substituted by fa + 1. 
Starting from fl23| , and using (j2ip and (T25l) . one finds: 



1 n 

N n ,n-{h+l) = — "2 ^ ^ njfc y4 A 



A /-^t ~*>n,k^k,n-{h+l) 

n ' n k = n-(h+l)+l 



1 - 



n-(h+l) 



1 1 - 

I ^n,n-^-n,n— (h+1) ~\ ^ ^ -^n,n— k-^-n— k.n— (h+1 



a h+ i 1 v-^ / 1 \- T-f ~ a 




, S IT ~ a h+l-k- 

" 1 peP(fe) n^ep u y 

For any partition p = {ri\,n2, ■ ■ ■ ,^( P )} £ 'P(^)) consider the set obtained 
by appending to the l(p) numbers rij in p the number n = (fa + 1) — k, 
that appears in the last factor (so that it can be rewritten in the form a no ). 
The set p' = {nQ,ni,ri2, ■ ■ ■ ,ni( P )} is now a partition of (fa + 1) of length 
l(p') = l(p) + 1. Considering then all possible values of k resulting from 
the summation of the last expression, we obtain in this way all possible 
partitions of (fa + 1) of length greater than 1. Considering the term outside 
the summation symbol as corresponding to the unique partition {fa + 1} of 
length 1 of (fa + 1), we can rewrite the last expression in the following form: 

N n ,n-(h+i) = — 52 ! [ 1 ■ 

a Q - LJ - a Q 

This is Formula ( 125]) written for the integer (fa+1) instead of fa, so our proof by 
induction is now complete and (125]) is true for all values offa = l,2,...,n — 1. 
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Usually it is more convenient to use the compositions in place of the 
partitions. Using the compositions, Formula (|25|) becomes 



N n , n . h = — Y] //(c) TT Vh = l,2,...,n-1, 

c€C{h) n 3 &c 

where the factor //(c) is given in Equation (J5]). 

We are now ready to calculate the elements v n from Equation 

n 

Vn = ^ Nn ' h Wh = 
h=l 

= I — Yl Vic) I f I (-a&ao) + — (-a n ao) 

h=l \ u ceC(n-h) 6 c u y u 

For any fixed h < n and c G C(n — ft.), the 1(c) + 1 numbers in c U {h}, 
obtained by appending the number h to the 1(c) numbers rij G c and reorder- 
ing, form a composition d of the number n of length 1(d) = 1(c) + 1. We 
can consider the last term in the above equation obtained by itself from a 
partition of the integer n of length 1, that is for the value h = n of the sum 
index. We can therefore write the previous formula in a more compact way: 

«»= e n — . 

d e C(n) rij 6 c' 

and this, with d substituted by c, is Formula ( ITU]) of Theorem [TJ Formula 
([9]) could be obtained in a similar way. 



3 Proof of Theorems El Eft SI and [5 



Proof of Theorem [2j We only prove Eq. (ITTj) concerning partitions. Eq. 
(j!2p has an identical proof. Simplifying the second member of Eq. (ITTjl . we 
obtain: 



"o e n 

p&V(n) rij ep 



a, 



«0 



-En n 



p&V(n) \rij£p 



rij gp 



a; 



a 
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e n(-^ )=(-D*» 



I! • 



p&V(n) nj£p 

where we have used Eq. that, by hypothesis, is true. The last member is 
equal to the first member of Eq. ffTTl) . so the theorem is proved. □ 



Proof of Theorem [3j Let n be any fixed positive integer. Using the formula 
for the binomial expansion one obtains: 



n 

' n 



o = (i-ir = E I i) in ~ h (- l ) h i 

that can be written in the following way: 

£(^^ = °- < 26 > 

Defining 

{-If 1 
ah hi ' an ~ h (n-h)\' 

Eq. (j2BJ) become like Eq. (J7J), with a Q = a = 1. Then Theorem [TJ applies 
and one may write: 

which is Eq. f[T3j) . The two other formulas in Theorem [3] follow immediately 
from Eq. (fl3|) taking out of the sum the contribution of the composition of 
length 1 of n and simplifying it with the first member of the equation. □ 



Proof of Theorem HI Let n be any fixed positive integer. We first prove 
Eq. ( 116]) . The first expression in Eq. (JTJ can be written in the following way: 

V - — = (27) 

^(n+l-h)\ h\ ' 1 ' 

Defining 



B, 



nl' n ~ n (n + l-h)\ 
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Eq. ( 127)) becomes like Eq. (J7J), with a = ct = 1. Then Theorem [T] applies 
and one may write: 

§ = £ Mc) n (-(^Tiyy 

which is Eq. (ITS)) . 

It is well-known that B 2n+1 = 0, Vn > 0. A purely algebraic proof of this fact 
is outlined in [TJ p. 238. It is based on a formula published by J. Faulhaber 
in 1631 and proved by B. Pascal in 1654, as reported in [6], Theorem 3.1 (i). 
Taking into account that B 2n+ i = 0, Vn > 0, we first prove that the Bernoulli 
numbers of even index B 2n satisfy the following relation: 

Y> 1 (2ft - l)B 2h _ 

^(2n + 2-2ft)! (2ft)! "°- (28) 

By hypothesis, all Bernoulli numbers of odd index greater than 1 vanish, so 
it is convenient in Eq. (127)) to separate the sums of the odd and of the even 
indices. We can write Eq. ( 127)) for n = 2m and for n = 2m + 1. When 
n = 2m, Eq. (127)) gives: 



2m R 2m 1 R R 
1 -Oft _ \ ^ 1 -Dft -Pi 

o (2m + l-ft)! ft! ~ ^ (2m + l-ft)! ft! (2m)! 

ft even 

1 B 2 k B\ 

(Om _L_ 1 _ 9.M' /^.M I 



, (2m + l-2ife)! (2Jfe)! (2m)! 



Si ^ 2m + 2-2A; B 



that is: 



When n = 2m + 1, Eq. (1271) gives: 



Ezm -|- Z — Zft. -D2A: 
(9™ 4- 9 - 9MI C9M! ^ ^ 



(2m)! ^ (2m + 2 - 2Jfe)! (2Jfe)! 



2m+l R 2m 1 R R 

El JDh \ 1 JDh JD\ 
= v + 
(Om 4- 9 - Ml 7)1 ^ (9<m4-9 — Ml 7)1 



( (2m + 2-ft)! ft! f-< o (2m + 2-ft)! ft! (2m + 1)! 

h even 



1 B 2k B x 

(Om 4- 9 - 97^1 (0h\\ 



(2m + 2-2k)\ (2k)\ (2m + 1)! ' 

E2m + 1 B 2 k 
(Om 4- 9 - OhW (OhW ' ' 



fc=0 

that is: 

Bx 2m + 1 B 



(2m)! ^ (2m + 2 - 2k) ! (2Jfe)! 
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The difference between Eqs. ( 15Uj) and (J2HJ gives: 

2m + 1 2m + 2-2k \ B 2k 



(2m + 2-2£;)! (2m + 2 - 2k)\ J (2k)\ 

m 

= E 



fc=0 

2fc - 1 ff 2fc 
, o (2m + 2 - 2Jfe)l (2fc)! 

Changing k with /i, this equation is identicai to Eq. ()28 
To prove Eq. ( ITT)) , one may define 

(2n - l)£ 2n 2 



(2n)! ' (2n + 2-2/i)! ' 

where the factors —1 in a n and 2 in a n -h have been introduced to have 
do = ozq — 1. Then Eq. ()2"g)) becomes like Eq. (jTj), with a = a = 1, and 
Theorem [TJ applies. By substituting a n and a n _^ just defined in Eq. flJDJ), 
one finds Eq. (ITTj) and the proof is complete. □ 



The proof of Eq. (1171) in Theorem H] is based on Eq. (1281) . This is an 
interesting result by itself, so we state the following Corollary, proved above. 

Corollary 1 The Bernoulli numbers of even index B2 n satisfy the following 
relation: 

A 1 (2h - l)B 2h = 

f^ o (2n + 2-2h)\ (2h)\ 



Proof of Theorem O Let n be any fixed positive integer. We first prove 
Eq. ( |18p . The first expression in Eq. (J2J) can be written in the following way: 



Defining 

_ En l + (-l)"- fe 

tth ~ n\ ' an ' h ~ 2{n-h)\ ' 

Eq. ( )3Tp becomes like Eq. (JTJ), with ao = olq = 1. Then Theorem [TJ applies 
and one may write: 



n 



f=E,(cn(-^F) 



ceC(n) nj&c 
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which is Eq. (fIB|). 

We next prove Eq. (1T9")) . It is well-known that the Euler numbers of 
odd index vanish. This can be seen also from Eq. (TJJ), in fact, in the sum, 
the factor [l + (— l) n_?i ] vanishes when n and h have different parity, and, 
starting from E\ = 0, obtained from (jlj), it is easy to realize that £>2n+i = 0, 
Vn > 0. Eq. (J2J) can then be written in the following way: 

n /o»A 

Eq = 1, E 2n+1 = 0, Yl[ 2 n h ) E2h = - ( 32 ) 



The third expression in Eq. ( 132)1 is also equal to 

1 E-jh 

(2n-2h)\ (2h)\ 



h=0 

Defining 



E2n 1 

a n -h 



n (2n)!' n ' n (2n-2/i)!' 

Eq. ( 133)1 becomes like Eq. (JJJ), with ao = ckq = 1. Then Theorem [T] applies 
and one may write: 

E^n \ ^ , , tt / I 



( 2; ,)! S II V (2//;i ! 

which is Eq. (jI5j). □ 



4 Conclusions 

Using only elementary results of linear algebra and some combinatorics we 
found the results expressed by Theorems [TH2l They allow to obtain explicit 
formulas for a wide class of sequences of numbers defined by recursion in 
terms of partitions of integers. We applied Theorem [T] to find formulas for 
the inverse of the factorials, the Bernoulli numbers and the Euler numbers. 
We think that the results expressed by Theorems [TH3] and by Eqs. (1X7)1 and 
(1T8)1 are new. Eqs. (1T61) and (Tl9)) have been also found by Vella [2] and 
Malenfant [3], but using tools from calculus. On the contrary our proofs are 
based only on some basic linear algebra and some combinatorics. 
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